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ABSTRACT 

The principal paradigm for the generation of the non-thermal particles that are responsible for the 
prompt emission of gamma-ray bursts invokes diffusive shock acceleration at shocks internal to the dy- 
namic ultrarelativistic outflow. This paper explores expectations for burst emission spectra arising from 
shock acceleration theory in the limit of particles cooling much slower than their acceleration. Para- 
metric fits to burst spectra obtained by the Compton Gamma-Ray Observatory (CGRO) are explored 
for the cases of the synchrotron, inverse Compton and synchrotron self-Compton radiation mechanisms, 
using a linear combination of thermal and non-thermal electron populations. These fits demand that the 
preponderance of electrons that are responsible for the prompt emission reside in an intrinsically non- 
thermal population, strongly contrasting particle distributions obtained from acceleration simulations. 
This implies a potential conflict for acceleration scenarios where the non-thermal electrons are drawn 
directly from a thermal gas, unless radiative efficiencies only become significant at highly superthermal 
energies. It is also found that the CGRO data preclude effective spectroscopic discrimination between 
the synchrotron and inverse Compton mechanisms. This situation may be resolved with future missions 
probing gamma-ray bursts, namely Swift and GLAST. However, the synchrotron self-Compton (SSC) 
spectrum is characteristically too broad near the vF v peak to viably account for bursts such as GRB 
910601, GRB 910814 and GRB 990123. It is concluded that the SSC mechanism may be generally incom- 
patible with differential burst spectra steeper than around E~ 25 above the peak, unless the synchrotron 
component is strongly self-absorbed. 

Subject headings: gamma-rays: bursts — radiation mechanisms: non-thermal — gamma rays: theory 
— relativity 
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are discussed in Section l2"lfl The focus is on GRBs observed by 
the Compton Gamma-Ray Observatory, specifically those with 
EGRET detections that provide broad spectral coverage and 
therefore more constrained fits around the spectral peak. It 
is found that acceptable fits are only possible with a marked 
dominance of non-thermal electrons, contrasting the particle 
distributions obtained in acceleration simulations. This poses 
a general problem for any burst acceleration paradigm where 
the non-thermal population is drawn probabilistically from a 
thermal gas. 

In addition, fits are found to be more or less equally accept- 
able for the synchrotron and inverse Compton processes; it is 
anticipated that Swift will provide a number of bright bursts 
where spectral fitting down into the X-ray band might enable 
discrimination between such emission mechanisms. However, at 
this juncture, spectroscopic discrimination against unabsorbed 
synchrotron self-Compton scenarios seems probable, given that 
the characteristically broad SSC spectra in vF v space are dif- 
ficult to reconcile with CGRO data. In service of this anal- 
ysis, a compact formalism for SSC emissivities from thermal 
and truncated power-law electron distributions is present in the 
Appendix. It is anticipated that GLAST's broad spectral cov- 
erage from ~ 10 keV to ~ 100 GeV, supplied by the Gamma- 
Ray Burst Monitor in the BATSE band and the Large Area 
Telescope in and above the EGRET band, will render GLAST 
a powerful tool in constraining and interpreting burst spectral 
properties in the future. 

2. SYNCHROTRON EMISSION AND ITS INTERPRETATION 

This section explores the relationship of synchrotron spec- 
tra to the electron energy and angular distributions, inferring 
Lorentz factor distributions using data from the bright burst 
GRB 910503 as a representative case study. 

2.1. Synchrotron Formalism and Electron Distributions 

The formalism for synchrotron radiation in optically thin en- 
virons is standard and readily available in many textbooks, in- 
cluding the treatments of Bekefi (1966), Jackson (1975) and Ry- 
bicki & Lightman (1979). Throughout this paper, the standard 
convention for the labelling of the photon polarizations will be 
adopted, namely that || refers to the state with the photon's 
electric field vector parallel to the plane containing the magnetic 
field and the photon's momentum vector, while _L denotes the 
photon's electric field vector being normal to this plane. 

The classical (photon) angle-integrated emissivities for the _L 
and || polarization for monoenergetic electrons of Lorentz fac- 
tor 7 moving at angle 9 to the field are (e.g. Westfold 1959; 
Jackson 1975; Rybicki and Lightman 1979) 
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in units of sec -1 , where the K v are modified Bessel functions 
of the second kind. The generally higher probability of emit- 
ting _L photons than || ones reflects the intrinsic dipolar radia- 
tive nature of accelerating/oscillating charges; the polarization 
is elliptical for photon propagation along the field. The char- 
acteristic synchrotron rate factor that scales these emissivities 
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Here, at = e /(he) is the fine structure constant. The photon 
energies e c that appear in these functions that form the scal- 
ing for the energies of emission are [with the electron Compton 
wavelength over 2ir being A = fi/(m e c)]: 
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where r L is the particle's Larmor radius, and B CI = 4.41 x 
10 13 Gauss is the quantum critical field. Note that since 7>1, 
the photon angles to the field form a narrow distribution about 
9 . Furthermore, the expressions in Eq. J5j are appropriate for 
cases where the particle momenta are not closely aligned with 
the magnetic field; situations of such near alignment are dis- 
cussed in Section 1231 below. 

The total synchrotron emissivity is obtained by integrating 
over the energy and angular distribution of the particles. For 
electron distributions n e (j, 9) in units of cm~ 3 , summing over 
polarizations forms the expression for the emissivity that will 
be used mostly throughout this paper (in units of cm~ 3 sec^ 1 ): 
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This is the differential photon spectrum; frequently in this work, 
reference will be made to the so-called vF v spectrum, a spectral 
energy distribution that is proportional to e 2 n s (e) . 

Physical diagnostics on the GRB environment are obtained 
by comparing the emission spectrum in Eq. Q with observa- 
tions. The obvious database of choice comprises BATSE obser- 
vations from the Compton Gamma- Ray Observatory (CGRO), 
specifically the BATSE spectroscopy catalog of Preece et al. 
(2000) with its pubicly-accessible electronic compilation. Yet 
this data possesses a principal limitation, namely that imposed 
by the restricted energy band of BATSE. Inferences of particle 
distributions are generally on a more secure footing when us- 
ing broader burst spectra, particularly in domains of significant 
spectra structure, which is the intrinsic nature of the BATSE 
pool. Thus, the two decades of energy afforded by BATSE of- 
ten are insufficient to describe the nature of the non-thermal 
particle distribution invoked in models. Data from the EGRET 
experiment on CGRO suitably augment this situation and ex- 
tend the dynamic range by 2-4 orders of magnitude in energy. 
Accordingly, a focus on EGRET bursts is insightful and of- 
fers more secure determinations of particle distributions. Use 
is therefore made of spectral compilations of BATSE, COMP- 
TEL and EGRET data that are offered in Schaefer et al. (1998). 
Specifically, a case study is made here of GRB 910503, observing 
that the inferences to be made generally extend to the handful 
of other EGRET bursts, and also to a large number of bright 
BATSE bursts that possess no confirmed EGRET detections. 

Note that the focus here is primarily on time-integrated spec- 
tra. Given the complexity and diversity of burst time profiles, 
and associated variations in spectral shape and the well-known 
hard-to-soft evolution, detailed spectral fits will depend on the 
time window chosen. A classic example for this is provided by 
the gradual emergence of a hard gamma-ray component in the 
spectrum of GRB 941017 (Gonzalez, et al. 2003). The main 
conclusions of this paper apply to the mean properties of the 



emitting particle distributions, and are striking enough to indi- 
cate their applicability for most or nearly all of the duration of 
the bursts studied herein. 

A variety of particle distributions are possible in modeling 
burst spectra. Clearly, all should contain a non-thermal compo- 
nent, since the data are clearly non-thermal: purely isothermal 
models do not work (Pacynski 1986; Goodman 1986), though 
constructed convolutions of thermal spectra of different temper- 
atures are possible. The emphasis here is motivated by physical 
interpretation, and the leading contender for non-thermal parti- 
cle energization in bursts is diffusive acceleration at relativistic 
shocks in the burst outflow. Such acceleration is well studied in 
non-relativistic flows, and is widely believed to originate with 
thermal particles that are diffusively transported when interact- 
ing with field turbulence in the shock environs. Such transport 
effects many crossings of the shock layer (e.g. see Drury 1983; 
Blandford & Eichler 1987; Jones & Ellison 1991 for reviews), 
that lead to a friction term in the momentum evolution, i.e. 
acceleration. While other acceleration models such as recon- 
nection and coherent electrodynamic energization exist and are 
appropriate in a variety of astrophysical environments (e.g. the 
solar corona, pulsar magnetospheres), diffusive shock accelera- 
tion is the most widely invoked in astrophysical models; it forms 
the centerpiece of discussion here as well as for many investiga- 
tions of bursts. 

Tavani (1996a, b) proposed a shock emission model that in- 
voked such diffusive acceleration as the means for generating 
non-thermal distributions. This work has been an interpreta- 
tive driver in the field, in no small part due to the impressive 
model fits of data obtained. The principal purpose of this Sec- 
tion is to delve deeper into the physical implications of such 
data fitting. Accordingly, it is appropriate to use an electron 
distribution similar to Tavani's ( 1996a, b) to add insight to the 
discussion. Here, the quasi-isotropic electron distribution 
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is adopted, where O(x) is a step function with O(x) = 1 for 
x > 1 and zero otherwise, and 7 T = kT/m e c 2 is the system or 
post-shock temperature. The factor e directly relates to the effi- 
ciency of acceleration, and is a parameter in the fitting protocol, 
as are 7 T (assumed much greater then unity) , the dimensionless 
77 , the non-thermal index S , and the normalization factor ng ■ 
The ensuing discussion will provide physical grounds justifying 
this choice of distribution in shocked environs; other possibilities 
with different mandates include the broken power-law adopted 
by Lloyd & Petrosian (2000). Note that ng is independent of 9 
for true isotropy. In many cases in this paper, ng ex 5(9 — ir/2) 
will be chosen for simplicity, since this corresponds to the dom- 
inant radiative contribution from the synchrotron mechanism. 

The parameter n is introduced to generalize from Ta- 
vani's (1996a, b) fits, permitting the non-thermal distribution 
to emerge perhaps at truly super-thermal energies, rather than 
exactly in the thermal peak (Tavani's restricted case). Typi- 
cal acceleration spectra in non-relativistic shocks assume forms 
similar in general appearance to Eq. (J5J. Specifically, they cor- 
respond to the n > 1 portion of parameter space (at least for 
ions), so that dissipation in the shock layer leads to a heating 
plus non-thermal acceleration that ensues only at super-thermal 
energies in the Maxwell-Boltzmann tail. Evidence for such a 
structured distribution can be found both in heliospheric mea- 
surements of accelerated ions near shocks (e.g. Earth's Bow 



Shock: see Ellison, Mobius & Paschmann 1990; interplanetary 
shocks: see Baring et al. 1997), and also a variety of acceleration 
simulations and theoretical analyses (e.g. see Scholer, Trattner 
& Kucharek 1992; Giacalone et al. 1993, Baring, Ellison & 
Jones 1993; Ellison, Baring & Jones 1996; Kang & Jones 1997). 
The situation for relativistic shocks is less certain, due in no 
small part to the absence of in situ spacecraft determinations of 
particle spectra; simulation expectations for such environments 
are discussed below. In this paper, values of n ~ 3 appear, 
being in general only weakly constrained by the fits and largely 
motivated by the above physical and theoretical considerations. 
While numerical evaluations of Eq. Q given the electron dis- 
tribution in Eq. JSJ were performed to obtain approximations 
to burst spectra, standard asymptotic limits are readily attain- 
able in analytic form. The non-thermal contribution to the syn- 
chrotron spectrum (from the truncated power-law), dominant 
whenever e ^> n 2+ , assumes the following power-law forms 
(derived, for example, from Eq. jllj): 
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Here, B±_ = B sin 9 is the component of the field perpendicular 
to the instantaneous electron motion. The thermal component, 
sampled when e <C n 2+ , necessarily declines exponentially at 
photon energies well above the thermal peak, but possesses the 
same power-law dependence as in the non-thermal case for low 
photon energies: 
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This spectral similarity at low energies follows from the narrow 
nature of the distributions being convolved with the synchrotron 
emissivity functions in Eq. J2J. Accordingly, it is not possi- 
ble to distinguish between emission from thermal or power-law 
distributions if spectrum measurements exist only at energies 
below the peak; some other diagnostic is required. The expo- 
nential decline above the thermal peak possesses a somewhat 
complex two-domain structure; the reader is referred to Pavlov 
& Golenetskii (1986), Brainerd & Petrosian (1987) and Baring 
(1988a) for the presentation of appropriate asymptotic forms in 
cases 7 T ^> 1 . 

It must be noted that throughout this paper, photon spectra 
are generated in the shock rest frame (SRF) and then boosted in 
energy to the observer's frame (OF) by a fixed Doppler factor, 
independent of photon angle in the SRF. This approximation, 
essentially adopted in many models of GRB spectra, is expedient 
but does not accommodate the fact that photon distributions in 
the SRF are inherently anisotropic due both to the instrinsic na- 
ture of the synchrotron and inverse Compton processes, and also 
to the fact that relativistic shocks naturally generate anisotropic 
ions and electrons. Notwithstanding, this boost approximation 
does not dramatically distort the angle-integrated OF spectral 
shape for typical SRF angular distributions and suffices to de- 
scribe the continuum for the purposes of the discussions here. 

2.2. Inferred Electron Distributions 

Before presentation of results, it is appropriate to address the 
modeling procedure used in this paper. Approximate spectra, as 
opposed to data fits, were derived visually using both log- log and 
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Fig. 1. — Photon spectra (a) from the synchrotron model and the corresponding normalized electron distributions (b) for the bright burst GRB 910503 
observed by the Compton Gamma- Ray Observatory. The data compilation presented in (a), adapted from Baring (1995), is a pre-publication release of 
that in Schaefer et al. (1998). The operating energy bands of the three detecting instruments, BATSE, COMPTEL and EGRET, are as indicated. In 
panel (b), the distributions are normalized to unit area using a scaling of the Lorentz factor in terms of 7 T (see text). The solid curves correspond to 
the particle distribution/resulting synchrotron spectrum pair that approximates the observed continuum very well [labelled synchrotron 'fit' in panel (a)]; 
model parameters included r\ = 3.0 and e ss 2550. The dashed case (for r\ = 3.0 and e = 48.4) illustrates a thermally-dominated situation that fails to 
account for the observations. In both cases, the non-thermal distribution index was 8 = 3.6 , the electron pitch angle was assumed to be 9 = 7r/2 , and 
T'y%,B±/Bcr = 0.90 arose from adjustments along the energy axis. The open triangles depict the approximate distribution obtained from the full PIC 
simulation of a plane- parallel relativistic electron-positron plasma shock (Mace and Jones 1994; see text). 



semi-log representations of spectra to increase confidence in the 
results. Least squares trials were performed in some cases, but 
with no addition of insight or perceived accuracy. This choice 
of a visual approach was one of expediency, and did not com- 
promise the scientific conclusions presented here; the essential 
points are unambiguously made without the need of statistical 
analysis. This latitude is afforded by the limited energy range 
of the data in the BATSE catalog: slight adjustments in param- 
eters used to fit the extant data for a given burst would almost 
certainly have been dominated by those needed to accommodate 
a broader spectral range. For the specific cases of GRB 910503 
and GRB 910601, where EGRET data are incorporated in the 
study, the uncertainty in the > 100 MeV data points is a major 
contributor to the determination of model parameters. In bursts 
where no EGRET or COMPTEL data were used, the "fit" was 
actually to the Band (1993) model spectral fit presented in the 
Preece et al. (2000) spectroscopy catalog, not to the actual data. 
In such cases, the high energy Band model index was adopted, 
and other distribution parameters adjusted. In this sense, there 
was in implicit second order incorporation of the statistics em- 
bodied in the Band model fit. It must be emphasized that there 
is no uniquely preferred method for statistically fitting public 
data since it is already a convolution of model and instrumental 
response. The correct approach (as in Lloyd & Petrosian 2000) 
would be to fold spectral models here through the BATSE or 
EGRET detector response function, an enterprise that rapidly 
becomes more ambitious when different instrumental responses 
must be blended. Such complexity seems premature for the ex- 
plorations of this paper, and is deferred until the more refined 
data of Swift and GLAST becomes available. 

Results for the spectral approximation of the GRB 910503 
data are presented in Figure ^ together with the multi- 
instrument, time-integrated CGRO data as provided by Schaefer 



(private communication) prior to the public release in Schaefer 
et al. (1998). We remark that the spectral deconvolutions pub- 
lished in Schaefer et al. (1998) were made independently for 
each data set from the different CGRO instruments, applying 
an artificial normalization to reconcile offsets in the different 
data amplitudes. Such a procedure is less accurate than a joint 
spectral deconvolution of the entire spectral data set (see for ex- 
ample the analysis of GRB 990123 data in Briggs et al. 1999). 
Notwithstanding, the major conclusions of this paper are robust 
and are not affected by such subtleties of data analysis. 

Figure ^ displays both synchrotron emission spectra for the 
case of pitch angles 9 — 7r/2 (the so-called quasi-isotropic case), 
and the inferred electron distribution that generates the contin- 
uum. For the solid curve approximate case, the overall nor- 
malization is obviously a free parameter and was offset slightly 
from the data for visual clarity. The electron distributions were 
normalized to unit area in 7/7T space. The solid curve de- 
picts the case that would be required to explain the photon 
spectrum. The non-thermal particle distribution clearly dom- 
inates the thermal component in this "fit." Its index differs 
from the value of 3.8 that Tavani (1996a) obtained in his anal- 
ysis of this burst, a discrepancy that is actually insignificant 
given the error bars in the EGRET data. Adjustment of the 
rj parameter to improve the appearance of the model approxi- 
mation around and above the spectral peak led to the assumed 
value of 77 = 3.0 , accurate to around 15%. Notwithstanding, 
there appears to be some spectral structure around the peak 
that cannot be explained by the simple particle distribution in 
Eq. (J5J • The positioning of the peak in energy produced the fit 
value Tj^,B±/B cr = 0.90 , where T is the Lorentz factor asso- 
ciated with the bulk motion of the GRB expansion. Here, B is 
interpreted as the magnetic field in the burst expansion frame, 
i.e., not in the ISM or shock rest frames. 



The parameter e describing the relative normalization of the 
non-thermal and thermal components was found to be e « 2550 
( ± 10%), a strikingly large value that is not addressed explicitly 
in detail by Tavani (1996a). A discussion of the implications of 
an inferred particle distribution with such a large e follows just 
below. A similar property is obtained for synchrotron model 
fits for BATSE bursts GRB 940619 (trigger number 3035_6) 
and GRB 940817 (trigger number 3128_5) from the Preece et 
al. (2000) catalog, and also for GRB 910601 and GRB 910814 
whose composite spectra are presented in Schaefer et al. (1998). 
The spectra for bursts GRB 930131 and GRB 940217 that pos- 
sess EGRET detections also suggest a similar result, though 
the lack of a composite broad-band multi-instrument spectrum 
for each in the literature precluded a detailed analysis. The 
parameter regime of e > 1 is expected to be pervasive for 
BATSE bursts, due to the smoothness of their circumpeak con- 
tinua. For comparative purposes, another distribution and its 
resultant synchrotron is illustrated in Figure^as dashed curves. 
The distribution was deliberately chosen to possess a somewhat 
dominant thermal component, motivated by expectations from 
shock acceleration theory, and clearly yields a spectrum that 
grossly deviates from the observations both well above and well 
below the peak. While these deductions are obtained using a 
time-integrated spectrum, such a conclusion must apply either 
for most of the burst duration, or at sub-intervals when the 
burst was at its brightest. 

For the synchrotron mechanism to be a viable explanation of 
GRB spectra, the underlying particle distribution must be phys- 
ically realistic. Here, expectations of particle distributions from 
shock acceleration theory are the predominant focus, partly be- 
cause this was the preferred acceleration mechanism in Tavani 
(1996a, b), and largely because this is the underlying energiza- 
tion assumed in most burst models. In a nutshell, both theory 
and observation pertaining to non-relativistic shocks strongly 
indicate 77 > 1 and e >• 100 cases. Note that e/77 4 is essen- 
tially the parameter that demarcates the relative normalization 
of the non-thermal and thermal components in Equation (jSJ) . 

Applications of shock acceleration in the heliosphere and at 
supernova remnant shells more or less mandate these parameter 
regimes, since the seed pool of particles that can be subjected 
to acceleration are thermal and moreover plentiful. In hclio- 
spheric environs, the source is the cold solar wind, and there is 
little evidence of significant pre-heating that is not connected to 
the shock system under study. The in situ spacecraft measure- 
ments exhibit strong heating at shocks and the generation of 
dominant thermal ions in the downstream regions, correspond- 
ing to e/77 4 < 1 cases. Data taken at the high sonic Mach num- 
ber Earth's bow shock (e.g. see Ellison, Mobius & Paschmann 
1990; Scholer, Trattner & Kucharek 1992) and low sonic and 
Afvcnic Mach number interplanetary shocks (e.g. see Gosling 
et al. 1981; Baring et al. 1997) clearly support this assertion. 
In supernova remnants, direct particle measurements are, or 
course, impossible, and inferences can only be made from ra- 
diative signatures. While electron temperatures and densities 
can be estimated from X-ray signals and molecular emission in 
the IR/optical bands, the data cannot yet clearly discriminate 
e/77 4 <§; 1 and e/77 4 ^> 1 situations. 

From a theoretical perspective, diffusive shock acceleration 
naturally draws particles from the thermal pool via transport 
in shock-associated turbulence. Whether this be first or second 
order in momentum transport, thermal and non-thermal par- 
ticles are subjected to similar diffusion, at least as far as ions 
are concerned, so that e/77 4 > 1 circumstances would require 



virtually complete suppression of thermalization and a 100% 
efficiency of acceleration out of the thermal pool. Such a sce- 
nario would require unusually anomalous spatial diffusion co- 
efficients k(p) (i.e. extremely small) at thermal momenta, or 
perhaps strange and non-dissipative electrodynamic properties 
in the shock layer. Note that n(p) is more likely to increase 
at low energies due to the possibility of wave damping by ther- 
mal populations (e.g. see Forman, Jokipii & Owens 1974, for 
the relationship between k(j>) and wave power spectra within 
the confines of quasi- linear theory). Without the assumption 
of anomalous forms for n(p) , Monte Carlo transport simula- 
tions (e.g. see Baring, Ellison & Jones 1993; Ellison, Baring 
& Jones 1996) clearly indicate a downstream thermal popula- 
tion that dominates the extended non-thermal distribution in 
number, corresponding to e/77 4 < 1 ■ The "injection" of ions 
from the thermal pool is extremely efficient, but not incredi- 
bly so. Hybrid plasma simulations (e.g. see Trattner & Scholer 

1991, 1993; Scholer, Trattner & Kucharek 1992; Giacalone et al. 

1992, 1993; Liewer, Goldstein & Omidi 1993) of wave generation 
and ion transport near non-relativistic shocks, which treat the 
electron component as a background fluid and therefore address 
only some of the electrodynamic aspects of the plasma in the 
shock layer, produce similar e/77 4 < 1 distributions. 

These indications are strongly suggestive, but not decisive 
for the gamma-ray burst problem, since they largely don't fo- 
cus on electron properties and equally importantly do not ex- 
plore acceleration at relativistic shocks. In situ measurements 
of non-thermal electrons at shocks are limited; examples include 
recent Geotail observations of interplanetary shocks (Shimada, 
et al. 1999; Terasawa et al., 2001). Accordingly, discussion is 
dependent on simulational information, as is the case for rel- 
ativistic shocks, and the literature is sparse on such a subject. 
Monte Carlo simulations of ion acceleration at relativistic shocks 
(e.g. Ellison, Jones & Reynolds 1990; Baring 1999; Ellison & 
Double 2002) generate the same injection properties as do non- 
relativistic shocks, i.e. e/77 4 < 1 • 

Plasma simulations of relativistic shocks are similarly sparsely 
studied. Gallant et al. (1992) and Hoshino et al. (1992) pre- 
sented results from one-dimensional (ID) particle-in-cell (PIC) 
simulations of electron-positron and electron-positron-proton 
plasmas in ultrarelativistic flows with perpendicular shocks, 
where the mean field was orthogonal to the shock normal. Such 
simulations treat the full electrodynamics of such plasmas, but 
are severely limited by CPU constraints. Accordingly, they tend 
to use unrealistically small proton to electron mass ratios, and 
model relatively small spatial scales. These constraints tend 
to mask the production of some wave modes and limit consid- 
eration to particles with energies at most a few to ten times 
thermal. The context for these investigations was pulsar wind 
flows, and they revealed no acceleration in pure e plasmas. 
With the addition of protons in significant abundances, Hoshino 
et al. (1992) demonstrated that non-thermal positrons (and not 
electrons) could be created, as the ions generated elliptically- 
polarized magnetosonic modes that resonantly interacted with 
the positrons. 

Mace and Jones (1994, unpublished; F. Jones, private commu- 
nication) explored a ID PIC simulation similar to Gallant et al. 
(1992) for an ultrarelativistic e ± plasma shock of arbitrary field 
obliquity, and found that electron/positron acceleration arose in 
quasi-parallel scenarios, but disappeared when the field obliq- 
uity to the shock normal increased above around 0bn ~ 30 
degrees, concurring with zero acceleration in the earlier results 
for perpendicular shocks. A scaled version of their parallel shock 



e distribution (binned in energy) is depicted in Fig.^ and is a < —2/3 clearly marked. The most common fitting model 

used in the catalog was the popular Band (1993) GRB model 
for the observed flux / 7 (e) in photons/cm 2 /sec: 



fairly proximate to the r\ = 3 , e = 48.4 case. A notable feature 
of the simulation data is the very limited range of energies ac- 
cessible to the PIC technique, imposed by the severe CPU lim- 
itations for such complex codes. Some interesting recent PIC 
simulations results have been presented by Shimada & Hoshino 
(2000) and Hoshino & Shimada (2002) for electron-proton plas- 
mas containing mildly-relativistic shocks, where suprathermal 
power-law electrons are generated and the distribution closely 
resembles the PIC simulation data depicted in Fig. ^ Yet, as 
with the Mace and Jones results, the complete distribution re- 
veals that the power-law emerges only well into the Maxwellian 
tail, corresponding to the e/r] 2+s <J 0.1 situation elicited in the 
Monte Carlo simulations. 

It is noted in passing that the one-dimensional nature of these 
simulations potentially imposes an important absence of critical 
physics to the acceleration problem: Jokipii, Kota & Giacalonc 
(1993) and Jones, Jokipii and Baring (1998) demonstrate that 
particle diffusion across field lines is eliminated when simula- 
tions restrict the spatial dimensions to less than three. The 
implication of this result for the e ± shocks is a cessation of ac- 
celeration when thermal particles cannot be transported across 
fields to return to the shock from the downstream side; Bar- 
ing, Ellison & Jones (1993) demonstrated that for high Mach 
number, non-relativistic ion shocks, a quenching of acceleration 
through this effect would arise at around 9bn ~ 30° , a result 
that matched the behavior seen in the rclativistic PIC simula- 
tions of Mace and Jones. Moreover this result is consistent with 
the absence of acceleration in non-relativistic ID hybrid plasma 
simulations applied to the quasi-perpendicular hcliospheric ter- 
mination shock (e.g. Liewer, Goldstein & Omidi 1993; Kucharek 
& Scholer 1995), a site that is believed to generate the non- 
thermal anomalous cosmic rays. The essential loss of physics in 
simulations of restricted dimensionality is clearly crucial to the 
correct modeling and interpretation of particle acceleration in 
gamma-ray bursts. 

The principal conclusion of this discussion is that there is little 
observational or simulational evidence for quasi-isotropic elec- 
tron distributions in the environs of shocks that would render a 
synchrotron spectrum commensurate with the GRB910503 data 
in particular and typical burst spectra in general. 

2.3. The Low Energy Spectral Index Issue 

An issue that has already been substantially discussed in the 
literature is whether the synchrotron model is compatible with 
spectral slopes below the burst peak. For GRB 910503, Fig. ^ 
clearly indicates that there is no conflict between the data and 
the model. The fixed value of a = —2/3 for the low energy 
index in Eqs. and [7| suggests that this is not the situation 
for all bursts, a fact that was demonstrated by Preece et al. 
(1998; see also Preece et al. 2000; 2002). Synchrotron theory 
distinctly predicts this canonical —2/3 index at low energies 
for any isotropic electron distribution that has a characteris- 
tic Lorentz factor 7 C as its lower scale. Correlations between 
Lorentz factors and pitch angles can alter this property, as dis- 
cussed in Section 12 . 51 below. Any index greater than this value 
cannot be accommodated by quasi-isotropic distributions; in- 
dices less than this are easily described by a steepening particle 
distribution above j c . Preece et al.'s (2000) display of the dis- 
tribution of low energy indices clearly indicated a conflict, as 
illustrated in Fig. There the low energy indices for various 
empirical fitting models used in the BATSE spectroscopy cat- 
alog are binned in a histogram, with the permitted range of 
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Most bursts possess values a < and (3 < . In folding the 
Band model, broken power-law models and others through the 
BATSE detector response, Preece et al. (2000) found that the 
Band model performed best in roughly 2/3 of the burst sample. 
Note that the index a recorded in the catalog is the asymp- 
totic low energy index (Preece et al. 2000), and does not ex- 
actly match the index local to the lowest energy ~ 25 keV in 
the BATSE window. The inferred values of a are subject to in- 
strumental selection effects, particularly when the peak energy 
is low. Hence, there is the possibility that the a distribution 
might move to slightly higher a when the detection band moves 
to lower energy, the situation for Swift. Note that such an effect 
was observed by Preece et al. (1998) when applying a correction 
to accommodate the energy of the BATSE window. A depiction 
of the a histogram for time-integrated spectra in the BATSE 
spectroscopy catalog is given in Lloyd & Petrosian (2000); it 
exhibits a similar distribution, as expected. 
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Fig. 2. — The low energy power-law index a distribution presented in 
Figure 7 of the BATSE spectroscopy catalog of Preece et al. (2000). The 
index is obtained by Band model or broken power-law model fits to burst 
spectra (see text). The public database consists of some 5500 spectra from 
sub-intervals in 156 bursts, for which a subset of 4776 spectra possess well- 
determined a . Domains where a quasi-isotropic synchrotron (a < —2/3; 
75.9% of spectra) or inverse Compton ( a < ; 97.8% of spectra) models 
can account for the observed low-energy index are as indicated. 

The quasi-isotropic synchrotron model clearly cannot account 
for nearly 25% of sample burst spectra. This is true even 
for pathological choices of electron energy distributions: from 
the form in Eq. Q, it is easily shown that the spectral index 
a = — 91og 10 n s (£)/dlog 10 e is a monotonically increasing func- 
tion of energy e for any electron distribution, due to the broad 
and smooth nature of the synchrotron emissivities in Eq. J5|. 



While this is disturbing, it may be more of an indication of 
problems with assumptions underpinning the a < —2/3 bound 
than a suggestion that synchrotron radiation is not responsible 
for GRB prompt emission. For instance, relaxing an adherence 
to isotropy or ng oc 6(9 — 9 C ) distributions can potentially alter 
the spectral slope. It is demonstrated in Section |2~11 below that 
for an ng independent of 7 , the spectrum below the peak can- 
not be flattened beyond a = —2/3 cases, and in fact steepens. 
However, introducing a 7 -dependent angular distribution can 
suitably ease the a bound; this is exemplified in the focus on 
pitch angle synchrotron emission in Section 12.51 

An alternative recourse is synchrotron self-absorption, i.e. re- 
laxation of an optically thin assumption. Given the theoreti- 
cal expectation of extremely flat spectra below the r syn = 1 
break ( eh s ex e 5 ' 2 for power-law electrons and a eh s oc e 2 
Raleigh- Jeans tail for thermal e~ : see, e.g., Rybicki & Light- 
man 1979), different groups have invoked such a mechanism 
to account for a > —2/3 cases in the BATSE catalog (e.g. 
Crider & Liang 1999; Granot, Piran & Sari 2000). Pushing 
the absorption break to gamma-ray energies requires extreme 
physical conditions (n e ~ 10 8 cm~ 3 , B ~ 10 8 Gauss), as noted 
by Lloyd & Petrosian (2000), and can be inferred from Eq. ©. 
These parameter regimes would be more extreme in the absence 
of the strong Doppler boosting associated with GRB outflows. 
Equipartition is usually invoked to justify such high fields, how- 
ever a widely-accepted and well-understood mechanism for the 
generation of such fields remains elusive. A convenient resolu- 
tion of this issue is afforded by the proposition of Panaitescu 
& Meszaros (2000), namely that the prompt emission is an in- 
verse Compton image of a self-absorbed synchrotron spectrum. 
The concomitant moving of the absorption break energy e t to 
below the optical band lowers the required values of B and 
n e to less extraordinary ranges. The X-ray spectral indices in 
the scenario of Panaitescu & Meszaros (2000; see also Liang, 
Boettcher & Kocevski 2003) can comfortably accommodate the 
a = domain, as discussed in Section l3.4l 

An additional concern is that the turnover energy e t (in the 
GRB expansion frame) satisfying r syn (e t ) = 1 is a function 
of a number of environmental parameters. Equation (6.53) of 
Rybicki & Lightman (1979) for the optical depth r syn for self- 
absorption can be manipulated to arrive at the equation for St ■ 
This result can be found in many papers; here Eq. (7) of Baring 
(1988b) is adapted to accommodate the specific form for the 
non-thermal portion of the distribution in Eq. JSj : 
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Here, t t is the Thomson optical depth of the emission region, 
and B± = Bsm9, as before. Note that Y(x) represents the 
Gamma function in Eq. I|10|l . The density and size of the absorb- 
ing medium are incorporated in r T . These parameters possess 
implicit dependence on the bulk Lorentz factor of the flow, with 
different evolutionary dependences for adiabatic and radiative 
expansion cases (e.g., see Dermer, Chiang & Bottcher 1999). If 
Eq. © is interpreted in the frame of the burst outflow, then e t 
is in the EUV/soft X-ray band for typical burst Lorentz factors. 
Given the number of parameters (five in all) contributing to the 



computation of the self-absorption optical depth, it is difficult 
to comprehend how e t can be fine-tuned to the BATSE band in 
the observer's frame; such a conspiracy is an issue also for op- 
tically thin synchrotron emission, whose peak depends on three 
parameters, Y , 7 T and B±_ . 

A strong test of the self-absorption proposition could be of- 
fered by Swift, which should determine whether or not BATSE 
was fully capable of determining the "asymptotic" values of a , 
especially for those bursts with peaks around 100 keV or lower. 
The Gamma-Ray Burst Monitor (GBM) on GLAST will provide 
an additional broad-band probe. Should the BATSE a distri- 
bution be mirrored by the Swift and GBM ones, synchrotron 
self-absorption in the gamma-ray band could then only be op- 
crating in a small minority of bursts, rendering it less germane 
to the spectral shape discussion, except perhaps for synchrotron 
self-Compton models. 

2.4. Influences of Anisotropy 

The foregoing discussion has focused on the case of syn- 
chrotron emission from quasi-isotropic or large pitch angle elec- 
tron distributions. Distributing the angles of electrons in a non- 
uniform manner can alter the spectral index of the emitted ra- 
diation. A prime example of this was the flattening incurred 
in distributing magnetic pair creation turnovers in neutron star 
environments, probed by Baring (1990). Here, obviously a dis- 
tribution skewed towards small pitch angles is a potential can- 
didate for modifying a . Such a distribution is not at all con- 
trived: synchrotron cooling can produce it quite naturally (e.g. 
see Brainerd & Lamb 1987). However, as alluded to above, it is 
now demonstrated that such distributions, if uncorrelated with 
Lorentz factor, yield only spectral steepening, and so cannot 
resolve the low energy index problem. 

The spectral behavior can adequately be demonstrated us- 
ing the non-thermal portion of Eq. JSJ, which for the pur- 
poses of this discussion will be written in the form n e (pf, 9) = 

fg n e (6 - 1) 7^" 7~ d for 7 > 7 min and is zero otherwise. 
Here 7 m i n = ?77t ■ If £ B = B / B cr is defined as the dimension- 
less cyclotron energy, then the differential rate of production of 
synchrotron photons in Eq. J3J can be quickly manipulated by 
reversing the integrations over 7 and y . The result is 
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This form for the synchrotron rate will be used later in the dis- 
cussions on synchrotron self-Compton emission. By adding and 
subtracting suitable forms like Eq. I jllj l with different choices 
of 7min , a general expression for the synchrotron emission 
from a piece-wise continuous power-law distribution of elec- 
trons can be obtained. Defining the photon spectral index to 
be a = — 91og 10 7T-s(e)/(ilog 1 oe, it is easily shown that a is a 
monotonically increasing function of energy e for any electron 
distribution with a low 7 cutoff, using piecewise continuous 
power-laws to successively refine approximations to a smooth 
distribution. 



Eq. (|llfl possesses two readily identifiable asymptotic limits, 
namely when e <C £B7min sm ^ an ^ £ ^ e B7min sm ^- These 
lead to the forms given in Eq. (JHJ • Now consider an angular dis- 



tribution defined over the range 9 U 
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The superposition of spectral forms in Eq. 10 then immediately 
indicates that such an angular distribution must generate a spec- 
trum with asymptotic forms e~ 2 / 3 at e <C £ B 7mi n sin6> m i n , and 
£ -( , 5- 1 )/2 w } ien £ ^ £ B 7min sm ^max • The spectrum in between 
depends on how rapidly fg varies with angle, and can be either 
of slope 2/3 or steeper, but never flatter than e~ 2 / 3 and never 
steeper than e~( |5_1 )/ 2 . This follows from the convolution over 
angles of spectra that are monotonically declining in energy nec- 
essarily being also monotonically declining in e . This assertion 
holds even for fg that is strongly peaked towards 6* min , and 
monotonically declining with 9 . Hence, angular distributions 
with 9 uncorrelated with 7 cannot ameliorate the low energy 
index issue due to the spectral shape implied by Eq. (II IB - 

2.5. Pitch Angle Synchrotron 

A more viable possibility for flattening the low energy in- 
dex is when there is a correlation between the pitch angle 9 
and the Lorentz factor 7 of charged particles. This provides a 
distinct kinematic difference from the foregoing considerations 
of synchrotron radiation from quasi-isotropic electron distribu- 
tions. The archetypical example for consideration here is so- 
called pitch-angle synchrotron radiation, a mechanism explored 
in detail by Epstein (1973; for applications, see Epstein & Pet- 
rosian 1973). This mechanism corresponds to synchrotron emis- 
sion from particles possessing very small pitch angles 9 <C I/7 . 
The polarization-averaged emissivity derived by Epstein (1973) 
in this limit, for electrons in a uniform magnetic field, can be 
written in the form 
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This form is simply obtained from the flux emissivity in Equa- 
tion (10) of Epstein (1973), and is applicable for e <; e c /j ; at 
higher energies, the usual synchrotron formula in Eq. (0} is op- 
erable if 7$ ^> 1 . Note that this process exhibits strong circular 
polarization due to the near coalignment of particles with the 
field, though there is significant linear polarization in narrow 
bands of photon energy. Observe that ^f = 37s/ (4e c ) <; 1 de- 
fines the intrinsic energy scale of pitch angle synchrotron (PAS) 
emission. Eq. (|14|) is virtually independent of e for domains 
$< 1, thereby defining an extremely flat spectrum. This prop- 
erty is a direct consequence of the kinematic coupling 9 ~ I/7 , 
yielding a dependence e ~ jB/B cr that contrasts the isotropic 
synchrotron dependence e ~ -f 2 8B /B cr . This weaker depen- 
dence on 7 immediately implies an inference of stronger source 
fields in order to move the characteristic PAS photon energy 
into the BATSE window. 

The spectrum drops dramatically above \P ~ 1 if j0 -C 1 , 
otherwise it merges into the e -2 / 3 synchrotron form if 76* ;> 1 . 
These shapes imply (i) that prompt GRB spectra that are flatter 
than e° below the peak (i.e. around 98% of bursts: see Fig- 
ure^ can be easily accommodated by a pitch angle synchrotron 
scenario (just like the inverse Compton scattering case consid- 
ered in Section I5TT1 below') . and (ii) accordingly that a variety of 
GRB spectral shapes can be accurately modeled due to the char- 
acteristically narrow peaking of the PAS emissivity in Eq. I|14|) . 
These were conclusions of Lloyd- Ronning and Petrosian (2002), 
who explored fits to BATSE bursts for pitch angle synchrotron 



models, and addressed correlations between spectral parameters 
and temporal properties. Hence, pitch angle synchrotron is an 
attractively versatile spectroscopic possibility for bursts. Note 
that the mechanism of "jitter" radiation in turbulent magnetic 
fields (Medvedev 2000) possesses similar spectral properties due 
to the intrinsic 9 ~ I/7 kinematic coupling, thereby offering an 
alternative possibility of similar promise. 

Pitch angle synchrotron might arise when electrons stochas- 
tically increase their pitch angle from virtually zero during an 
acceleration process. For it to provide the major contribution, 
a prevalence of small pitch angles must be strongly favored, i.e. 
the pitch angles 9 ;> 1/7 must not be significantly populated, 
so that normal synchrotron emission does not dominate. It is 
not clear how this might come about, though rapid cooling and 
inefficient diffusion transverse to the field is an interesting pos- 
sibility (e.g. see Lloyd- Ronning and Petrosian 2002). In addi- 
tion, small pitch angles must be prepared in the first instance, 
which suggests a coherent electrodynamic mode of acceleration 
as opposed to diffusive one. While PAS may afford a means of 
enhancing the non-thermal radiative signal relative to the ther- 
mal component, a detailed exploration of PAS spectral shapes 
and polarization for thermal and non-thermal distributions in 
the context of diffusive acceleration will be the subject of future 
work. 



2.6. Cooling Issues 

While detailed exploration of cooling effects is beyond the 
scope of this paper, a brief discussion is appropriate. Cooling 
is known to strongly influence spectral shape in quasi-isotropic 
synchrotron scenarios (e.g. Ghisellini, Celotti & Lazzati 2000), 
introducing a second "line of death" with a value of a — —3/2 . 
This presents problems for models invoking strong cooling in 
generating spectra like those seen by CGRO. Such models are 
most applicable to GRB environments where the emission re- 
gions are somewhat remote from the acceleration ones. This 
is because even though synchrotron radiation is a remarkably 
efficient process, diffusive acceleration at shocks is even more 
so. Accordingly, coincident acceleration and emission zones of 
the same physical scale can exhibit prolific acceleration up to 
some maximum Lorentz factor 7 max , before radiative cooling 
in the magnetic field causes cessation of acceleration. The spec- 
tral turnover that marks this cessation is quasi-exponential in 
nature, and have not been observed in burst spectra to date. 
This combination acceleration/radiation scenario is widely en- 
visaged in the contexts of X-ray and gamma-ray emission from 
supernova remnants and blazars. Whenever multiple spatial or 
temporal scales are sampled, convolution of spectral contribu- 
tions can blur the turnover into a broad distribution, essentially 
the case explored by Ghisellini, Celotti & Lazzati (2000). 

The rate of acceleration scales as the inverse gyrofrequency, 
which has a weaker dependence on electron Lorentz factor 7 
than the synchrotron cooling rate. Specifically, the acceler- 
ation timescale (in seconds) can be approximated by r acc ~ 
0.1/ re ii?Tcv/(/3s n -B) , where Et c v is the energy of the electrons 
in units of TeV, /3 sn is the shock speed in units of c , B is the 
mean magnetic field strength (measured in Gauss) , and / ro i is a 
factor of order unity that accounts for corrections to the diffusive 
timescale appropriate for relativistic shocks. The synchrotron 
cooling timescale (also in seconds) is T syn ~ 300/(i? 2 £tcv) ■ 
Various versions of these formulae are widely available, though 
a specific discussion of them in the context of blazars is in Baring 
(2002); they are effectively applicable in the comoving frame of 



the burst outflow, as opposed to the observer's frame. The elec- 
tron energy and the magnetic field are coupled by placing the 
synchrotron vF u peak in the soft gamma-ray band, i.e. setting 
T 7 2 B - 10 14 Gauss. Here V ~ 10 2 - 10 3 is the bulk Lorentz 
factor of the outflow. These quickly lead to the ratio estimate 
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Acceleration is clearly much more rapid than synchrotron cool- 
ing for typical T > 10 2 , a property that changes only when 
considering electrons that emit above about 10 GeV, i.e. for 
r -f 2 B ;> 10 18 Gauss. In the special case of T ~ 1 , the result in 
Eq. 1)15(1 corresponds to the well-known synchrotron cooling cut- 
off property of the Crab pulsar wind nebula (e.g. see de Jager et 
al., 1996), namely that the spectral turnover occurs at around 
50 MeV, independent of the magnetic field strength B . 

3. INVERSE COMPTON IN THE GAMMA-RAY BAND 

Since significant questions surround the viability of the syn- 
chrotron mechanism in explaining burst spectra, consideration 
of alternative emission processes is warranted. Here, the focus 
is on inverse Compton scattering, another very efficient radia- 
tive mechanism. For the purposes of the discussion here, it 
comes in two relevant varieties: scattering off a narrow-band 
soft photon distribution, and upscattering of synchrotron pho- 
tons (synchrotron self-Compton, or SSC) that provide an in- 
herently broad seed spectrum. It shall be seen that the result- 
ing predictions for prompt emission spectra for these two cases 
are manifestly different. The inverse Compton (IC) mechanism 
possesses two obvious advantages over synchrotron radiation: 
(i) high values of B can be avoided, since there is no require- 
ment that the synchrotron peak appear in the gamma-ray band, 
and (ii) the low energy index can, in cases to be discerned in 
Sections 13.11 and 13.41 be as flat as a — . At the same time, 
inverse Compton scenarios (particularly the SSC case) have the 
disadvantage of requiring that synchrotron radiation or other 
components be unobservable at optical and lower wavebands 
during the prompt and early afterglow phases. The parameter 
constraints thereby imposed are discussed in Section 13.61 

The well-known criterion for IC to dominate synchrotron as 
an electron cooling mechanism (i.e. also in total bolomctric 
luminosity for a spatially uniform medium) is simply that the 
energy density u s of the soft photons exceeds that of the am- 
bient magnetic field, u B = B 2 /8ir (e.g. Rybicki & Lightman, 
1979). In general, the value of u s can vary enormously depend- 
ing on the particular invocation for soft photon generation. For 
the specific case of a SSC model, u s couples to the field and can 
be determined by integrating eh s (e) in Eq. Q over e and mul- 
tiplying by some fiducial photon escape timescale, t csc — R/c , 
that reflects the physical size R of the emission region. For 
monoenergetic electrons of Lorentz factor 7 m j n 3> 1 , the fa- 
miliar result (e.g. see p. 168 of Rybicki & Lightman, 1979) 
for the synchrotron energy loss rate P s can be used to write 
u s /u B = n e t csc P s /u B = 4cr T rf 0SC n e 7mm/ 3 • Hence, the SSC 
emission can dominate the bolometric luminosity for significant 
but not particularly large values of the product «. e 7mi n ■ ^ was 
precisely this property that led Gonzalez et al. (2003) to ar- 
gue that the intriguing hard gamma-ray component seen by the 
EGRET TASC in GRB 941017 was not due to SSC emission, 
but perhaps due to inverse Compton scattering seeded by an ex- 
ternal soft photon source. The criterion for an SSC component 
to fall in the soft gamma-ray band is r7 4 lin B / B CY ~ 1 , which 



can easily accommodate fields in the Gauss range for T ~ 300 
and 7„ lin ~ 100 — 10 3 . A more thorough exploration of SSC 
parameter spaces for burst prompt and afterglow emission can 
be found, for example, in Meszaros, Rees, & Papathanassiou 
(1994), Sari & Esin (2001) and Zhang & Meszaros (2001). 

3.1. Inverse Compton Formalism 

The formalism for the inverse Compton (IC) emissivity in 
optically thin environs is readily available in many astrophysi- 
cal textbooks, such as Rybicki & Lightman (1979). Since the 
description of IC polarization effects is more subtle than for syn- 
chrotron radiation, being dependent on anisotropics in the soft 
photon population, the developments here are restricted to cases 
where photon polarizations are summed over. The rate n ic (e) 
of emission, integrated over the angles of the upscattered pho- 
tons, is independent of the direction of momentum of the seed 
photons. Hence, no integration over 9 is required, contrasting 
Eq. Q for synchrotron radiation. However, the rate is depen- 
dent on the soft photon energy e s , which must be integrated 
over, yielding the form (in units of cto~ 3 sec^ 1 ) 
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djn e ( 1: 6) 4^ 



7 £s 



(16) 



Here a T is the Thomson cross section, no(e s ) is the energy 
distribution of soft photons, and 



f(z) = 2z\og e z + z + l-2z 2 @(z) 



4 7 2 £ S 



(17) 



for a Heaviside step function 9(z) = 1 when < z < 1 and 
8(z) = otherwise. 

Mirroring the treatment of synchrotron radiation, while nu- 
merical evaluations of Eq. (|16fl were performed using the elec- 
tron distribution in Eq. JSJ to obtain approximations to burst 
spectra in Section 13.31 asymptotic spectral indices are readily 
attainable. For monoenergetic soft photons, the non-thermal 
contribution to the inverse Compton spectrum (from the trun- 
cated power-law) is dominant whenever e ^$> rj 2+ , and assumes 
the following power-law forms: 



^ic(e) 



NT 



e 
7^£ 

v 7 2 £ S ; 



-(S+l)/2 



e s < e < 7^e s , 



7~e s < e. 



(18) 



The index at high energies, e ^> 7 2 £ s , is identical to that 
for synchrotron radiation, since both mechanisms emit pho- 
tons with the same dependence on the electron Lorentz factor: 
£ oc 7 2 e s • At low energies, the differential spectrum is approx- 
imately independent of the upscattered energy of the photon, 
due to the z — > limit of f(z) and the fact that electrons near 
the ?77 T cutoff dominate the low energy photon production. 
This index of a = has been known since the early days of 
application of inverse Compton to astrophysical problems: see, 
for example, Figure 3 of the comprehensive treatment of Jones 
(1968) that also illustrates that the Compton spectrum flattens 
further to e +1 in the domain e <J e s of downscattering. 

As in the case of synchrotron radiation, the thermal com- 
ponent, sampled when e ^C r j 2 +° ) necessarily declines expo- 
nentially at photon energies well above the thermal peak, but 
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possesses the same power-law dependence as in the non-thermal 
case for low photon energies: 



"10(e) 



TH 



ocr, e s < e < j*e s • 



(19) 



This follows because the low energy photons are produced by 
electrons mostly near the peak of the thermal distribution, again 
sampling f(z) « const. 

3.2. Synchrotron Self-Compton Formalism 

From the perspective of spectral shapes at and below 
the gamma-ray peak, it is important to explore the proper- 
ties of synchrotron self-Compton (SSC) radiation, to distin- 
guish its character from IC emission that is seeded by quasi- 
monoenergetic soft photons. Here the analysis will first be con- 
fined to the case of a purely non-thermal electron distribution, 
i.e. regimes e 3> rj 2+ , where the distribution can be written in 
the form n e (j, 6) = n e fg (5-1) 7^„ J~ 5 for 7 > 7 min and is 
zero otherwise. Hereafter, when making the identification with 
the non-thermal portion of the distribution function in Eq. J5J , 
7min = VIt > and one can set n e — > ne and fg — ► ej-j-r/ 1 ^ /(6— 1) 
to establish a correspondence. 

For such a truncated power-law, a compact analytic form 
for the SSC spectrum was obtained, essentially reducing the 
quadruple integration to one over a single variable. Defin- 
ing a fiducial photon escape timescale, t csc = R/c for a one- 
zone emission region of size R , the average density of syn- 
chrotron photons that seed the SSC emission can be written as 
^o(£s) = tesch s (£ s ) , so that the form in Eq. (|ll|l can be directly 
inserted into Eq. (|16|l . The analytic reduction is expounded in 
Appendix A, with the result being 



«ssc(e) 



= A s , 
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dyK 5/3 (y)Ps(y + / dyK 5/3 (y)Q s m 



(20) 



where 



iVssn = 



Ctf 



V3 



8,T 



A 7 



T T n e 

J — 



f 9 sine 



(21) 



is the effective scale for the rate of the synchrotron-self-Compton 
process. Here, r T = n e a T R is the Thomson optical depth of the 
emission region. The emergent photon energy is encapsulated 
in a single dimensionless parameter 



eB c 



6S±7n 



B± = Bsint 



(22) 



that denotes the energy scale of the peak for SSC emission, 
which putatively could be in the gamma-ray band, and is the 
only spectral structure in the emissivity. The functions Ps(z) 
and Qs(z) are relatively compact and are given in Eqs. i|A7() 
and l|A8(l in Appendix A in terms of elementary functions. 
Eq. I|2(J|I can be routinely integrated numerically. 

Asymptotic energy dependences of the SSC rate in Eq. (|20|l 
are readily obtained. In cases where x ^> 1 , the Qs term 
never contributes to the integration since it is exponentially sup- 
pressed by the behavior of the Bessel function. Accordingly, it 
is quickly deduced that the SSC spectrum approximately obeys 
h ssc (e) (x e^(' 5 + 1 )/ 2 at energies e ^> 6~/+ n Bj_/B cv , as expected: 



a convolution of an extended synchrotron spectrum with a trun- 
cated power-law electron distribution with index (6 + l)/2 is 
a power-law of index (5 + l)/2 . There is a slight logarithmic 
modulation as outlined in Eq. (|A12I) . 

The x ^ 1 limit can also be promptly obtained. The Ps 
term samples only the t/ <C 1 asymptotic limit of the Bessel 
function, i.e. the y -5 ' 3 form. From this, it is quickly verified 
that the y integration is proportional to x~ ■ It is simple 
to show using the Eqs. (|A7|) and (|A8|) in Appendix A that the 
Qs term possesses an identical dependence. The net result is 



(X £ 
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as embodied in Eq. (|A9|) . a property that has 



been evident in various SSC models of bursts (e.g. see Sari & 
Esin, 2001). These two limits can be summarized via 
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( 23 ) 
The low energy index images the synchrotron low energy tail, 

a result that is not surprising since such a tail is steeper (and 
therefore its constituent photons are more populous) than the 
tail resulting from IC upscattering of monoenergetic soft pho- 
tons in Eq. i(TB|) . Accordingly, the bound a < —2/3 is placed 
on SSC spectra so that this mechanism cannot aid with the low 
energy spectral index issues discussed in Section 12.31 

To facilitate the fitting of GRB spectra with SSC spectra, 
the case of ultra-relativistic thermal particles must also be con- 
sidered, for which the distribution can be written in the form 
n e (7, 8) = n e / e /(27 3 )7 2 exp{-7/7 T }. When making the 
identification with the thermal part of the distribution function 
in Eq. JSJ, one can set n e — > ng and fg — > 2^f T to establish a 
correspondence with correct normalization. In Appendix B, an 
expedient derivation of the emission spectrum is offered, making 
use of the result already presented in Eq. (|20|) for a truncated 
power-law electron distribution. Basically, the S — > 00 limit is 
taken to generate the spectrum for monoenergetic electrons, and 
this is then convolved with a relativistic Maxwcll-Boltzmann 
distribution. Setting 



N„ : 



Of 



V3 



16tt A 



-4-/0 sin 6» 



Tt 



(24) 



as the effective scale for the rate of the synchrotron-self- 
Compton process for the hot thermal particles, one arrives at 



M 



= A s: 



TH 



dyK 5/3 (y)J 



y 



1/4 



(25) 



where the function J(n) is defined in Eq. (|B3|I . and is steeply 
declining as k increases. Here, %t = £ B cr /[6B±j^\ . Asymp- 
totic forms for the spectrum in Eq. I|25|l are routinely derived, 
being presented in Eq. (|B5|I . As expected, an e -2 / 3 form arises 
at energies below the characteristic energy 6B±j^/B a , and the 
spectrum is a weakly declining exponential far above this. 

3.3. Electron Distributions for Inverse Compton 

Results for the spectral fitting of the GRB 910503 data for 
the case of inverse Compton scattering are presented in Fig- 
ure |3J The fits were performed using the formalism presented in 
Section 13.11 specifically with numerical evaluations of Eq. (|16fl . 
Again, the multi- instrument, time-integrated CGRO data as 
published in Schaefer et al. (1998) was used. The figure dis- 
plays both inverse Compton emission spectra for the case of an 
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CGRO Composite Spectrum for GRB910503 
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Fig. 3. — Photon spectra (a) from the pure inverse Compton model and the corresponding normalized electron distribution (b) for the bright burst 
GRB 910503 observed by the Compton Gamma-Ray Observatory. The data compilation and detecting instruments presented in (a) are as in Fig^ The 
solid curves correspond to the particle distribution/resulting inverse Compton spectrum pair that approximates the observed continuum very well [labelled 
Inverse Compton 'fit' in panel (a)]; model parameters included r\ = 2.63 and e ss 243. The photon spectrum was generated using monoenergetic soft 
photons of some energy eo . In (a) the thermal (labelled T) and non-thermal (labelled NT) components are illustrated by the dotted and dashed curves, 
respectively. In panel (b), the total distribution is normalized to unit area as in Fig0 The dashed curve in (b) is the non-thermal portion of the electron 
distribution, for which the index was <5 = 3.6 , and r7^.es = 1-32 arose from adjustments along the energy axis. 



isotropic, monoenergetic soft photon distribution, and the in- 
ferred electron distribution that generates the continuum. The 
solid curve depicts the case that would be required to explain 
the photon spectrum, with the overall normalization again be- 
ing a free parameter; this curve was offset slightly from the data 
for visual clarity. The electron distributions were normalized to 
unit area in 7/7T space. 

The non-thermal particle distribution clearly dominates the 
thermal component in this "fit," as is the case for synchrotron 
emission. The index 5 = 3.6 was controlled by the EGRET 
data, and so is identical to that in Fig. ^ Adjustment of 
the r\ parameter to improve the appearance of the model ap- 
proximation around and above the spectral peak led to the 
assumed value of n = 2.63 , again accurate to around 15%. 
The positioning of the peak in energy produced the fit value 
r^B±/B CI = 1.32 , where V is the Lorentz factor associated 
with the bulk motion of the GRB expansion. The parameter 
e describing the relative normalization of the non-thermal and 
thermal components was found to be e « 243 (±10%), still a 
large value though only a tenth of that obtained for the syn- 
chrotron fit. 

The conclusion is essentially identical to that in Section |2~2*1 
the observational data can only be accommodated by an electron 
distribution that has a preponderance of non-thermal electrons 
as Fig. 103). Such a determination was also obtained for BATSE 
bursts GRB 940619 (trigger number 3035.6) and GRB 940817 
(trigger number 3128_5) from the Preece et al. (2000) catalog. 
Moreover, the fits in Figs. ^ and |3| are of comparable quality, at 
least below 1 MeV and above 5 MeV; i.e. there is little potential 
in the CGRO data to discriminate between the two emission 
mechanisms. There is a small bump/excess in the theoretical 
spectrum in Fig. [31 just above 1 MeV that is a consequence of 
the sharper inverse Compton scattering kernel; its prominence 
can be muted somewhat by slightly reducing both e and the 
index S . Such a liberty may not be afforded in the GLAST 



era when burst spectral indices will be better constrained be- 
cause of its broad energy range, sampled by the Gamma-Ray 
Burst Monitor (GBM) in the BATSE band and the Large Area 
Telescope (LAT) in the EGRET band. 

The last point to note is that the inverse Compton fit is tol- 
erated here because the low energy asymptotic index is possibly 
not fully realized in the BATSE window. The data possess a 
slight curvature that may be consistent with a continued flat- 
tening of the spectrum below 30 keV, as would be expected 
for the inverse Compton mechanism. Yet the data match the 
synchrotron emissivity also, indicating that BATSE cannot dis- 
criminate between these two mechanisms on the basis of the 
soft gamma-ray portion of the bandwidth. This impasse should 
be addressed when Swift is launched in the coming year. Swift 
possesses the sensitivity in the critical 1 -100 keV band to probe 
the spectrum below the vF v peak and afford diagnostics on 
the emission mechanism. Furthermore, the spectral coverage 
afforded by the GBM and the LAT on GLAST will also provide 
powerful observational constraints. 

3.4. Low Energy Spectral Indices 

The low energy asymptotic forms in Eqs. (J18JI and (|19|l raise 
the possibility that inverse Compton emission can eliminate 
the so-called "line of death" problem that is pervasive in syn- 
chrotron models. The IC mechanism can accommodate low en- 
ergy asymptotic indices as flat as a — , with steeper indices 
easily being generated by appropriately distributed electrons. 
The domain of compatibility is indicated in Fig.[2J where around 
98% of BATSE burst spectra in the spectroscopy catalogue of 
Preece et al. (2000) can be accounted for by the inverse Comp- 
ton process. This is an enticing property, but is not a definitive 
statement of the operation of this mechanism in GRB prompt 
emission regions. It must be emphasized that the establishment 
of an a > domain is contingent upon the soft photons that 
seed the upscattering assuming a narrow energy (i.e. quasi- 
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monoenergetic) distribution. This caveat is exemplified by the 
SSC value of a being identical to that for synchrotron radia- 
tion due to the breadth of the synchrotron seed spectrum. Hence 
invoking a synchrotron self-Compton model only dramatically 
impacts the operating range of B and not the spectral shape 
in the very soft gamma-rays. 

Given the discussion in Scction l2. 41 pertaining to synchrotron 
radiation, it is natural to ask whether anisotropics in either 
the soft photon or electron populations can dramatically alter 
the spectrum. The inverse Compton photon production rate in 
Eq. (|16|l implicitly assumes isotropy of the soft photons consti- 
tuting no(e s ) . Furthermore, it assumes that the angular dis- 
tribution of electrons is not seriously deficient along the line of 
sight to the observer; i.e. that the viewing perspective is not 
towards the side of an electron jet. Calculating IC spectra for 
anisotropic distributions is quite involved, and moreover sensi- 
tive to the choice of geometry. Yet spectral behavior depends 
mostly on the interplay between such geometry and the kine- 
matics of IC scattering. 

If the electron distribution is highly collimated and the ob- 
server's direction lies outside the cone of jet collimation, then the 
kinematics of scattering requires a smaller angle # S catt of scat- 
tering in the electron rest frame to send the photon towards the 
observer. This small # S catt yields a lower final photon energy in 
the observer's frame, so that the emitted inverse Compton spec- 
trum is softer outside the jet core. The spectrum then flattens 
at lower energies (e.g. see Dermer & Schlickeiser 1993; Baring 
1994 for blazar contexts) due to the kinematic depletion of pho- 
tons in the observer's direction, and can be even slightly flatter 
than e° . A problem with this scenario immediately arises: the 
energetics are dominated by emission in the jet cone so that 
assuming "off-axis" IC emission places more severe constraints 
on the GRB radiative efficiency. Note that fluctuations in the 
jet orientation or structure would lead to the emergent inverse 
Compton spectra being highly variable both in flux (observed in 
bursts and therefore not a problem, in principle) and in spectral 
shape — any structure such as the vF v peak would be expected 
to vary in time rapidly in energy due to the scattering kine- 
matics. While modest fluctuations (as opposed to hard-to-soft 
evolution) are present in gamma-ray burst prompt emission, de- 
tailed modeling would be required to discern whether they could 
match those produced by anisotropic scattering scenarios. 

Any anisotropy in the soft photon population is not sampled 
until they become more collimated than the electrons; only then 
does the emissivity change appreciably from Eq. I|16[l . Then, 
electrons within the photon "cone" of beaming experience an 
almost full range of incident photon angles in their rest frames, 
so that they contribute an inverse Compton spectrum similar to 
that for isotropic soft photons. However, the majority of elec- 
trons lie outside this cone, and for these, collisions with photons 
are predominantly head-on, leading to a depiction of the softer 
photons from the overall inverse Compton spectrum just as with 
electron collimation. Hence, one expects that such soft photon 
anisotropics could in principal also flatten the spectrum below 
the vF v peak. 

3.5. Electron Distributions for SSC 

To provide a meaningful comparison with the preceding mech- 
anisms, GRB 910503 is again used as a focal point. Yet, so as 
to provide a variation, here a case study for another EGRET 
burst, GRB 910601 is also presented. Results for the spectral 
fitting for synchrotron sclf-Compton scattering are presented in 
Figure 0] The "fits" were performed using the formalism pre- 



sented in Section 13.21 specifically with numerical evaluations of 
Eqs. (|2"U|> and (|23|l . The figure displays the SSC emission spec- 
tra for the case of an isotropic electron distribution, for each 
burst. The solid curve depicts the case that would most closely 
explain the photon spectrum, with the overall normalization 
again being a free parameter. Due to the strong dependence of 
the characteristic energy on electron Lorentz factor, the values 
of T^B±/B cr required to move the SSC peak into the gamma- 
ray band lead to inferences of sub-Gauss fields if j T ~ T ;> 300 . 
Consider first GRB 910503. The most obvious features are 
that the "fit" is more difficult than for the emission mechanisms 
previously considered, with poor results for the lower energies 
much below the peak, and that the thermal contribution must 
be dramatically suppressed. The cause is easily identified: the 
SSC spectrum is intrinsically much broader than either the IC 
or synchrotron spectra. This is due to the characteristic energy 
for SSC embodied in the dimensionless parameter \ (see Equa- 
tion [52) scaling as 7^ (or Xt oc 7^ for the thermal case). 
Consequently, small changes in the electron 7 near the peak 
of the Maxwell-Boltzmann distribution incur large variations 
in the SSC photon energy so that the spectral peak is broad. 
This effect is manifested in the non-thermal SSC emissivity in 
Eq. H20fl . which leads to the excesses evident around 30 keV in 
the "fits" of Figure^] and more profoundly for the thermal SSC 
contribution, which is extremely broad. From Eq. (|B5() a notable 
conclusion can be drawn: that the thermal SSC spectrum in the 
depicted vF v representation only starts to fall above around 30 
MeV, and then only slowly, i.e. it is blueward of the peak from 
the non-thermal SSC contribution. This drives the requisite 
dramatic suppression of the normalization of a thermal com- 
ponent to the underlying electron distribution, leading to the 
same issues (but in this case more strikingly emphasized) for 
acceleration from thermal energies that were discussed above. 
Note that for values of e > 25 , a comparably satisfactory SSC 
fit can be obtained with S ~ 5 . Observe that for such a large 
S , the tail of the thermal Maxwellian emerges above that of the 
non-thermal distribution for a small range of Lorentz factors 

7 > 37t • 

The situation for GRB 910601 is even more constrained, due 
to the markedly steeper gamma-ray spectrum. The thermal 
electrons must be even less populous to accommodate the Comp- 
tel and EGRET data, and the non-thermal electrons must pos- 
sess a very steep index S = 8.0 , considerably steeper than the 
S = 5.5 required for a synchrotron fit. In fact, both the syn- 
chrotron and IC models had much greater success fitting the 
GRB 910601 data, though are not depicted. Steepening the in- 
dex to 6 ^ 10 produces only marginal improvement in the fit. 
The inevitable conclusion of this analysis is that the breadth 
of the SSC emissivity for even monoenergetic electrons (derived 
in Equation (|B6|I . and closely approximated by the S = 8.0 
case in Figure 0} is large enough to provide significant or severe 
problems for the SSC interpretation of prompt emission from 
gamma-ray bursts. This contention is not based solely on the 
GRB 910503 and in particular the GRB 910601 data; the simi- 
larity of the spectra of GRB 910814 (Schaefer et al. 1998) and 
GRB 990123 (Briggs et al. 1999) to that of GRB 910601 strongly 
suggests that a spectroscopic incompatibility of the SSC mecha- 
nism for prompt emission may be commonplace, specifically for 
differential burst spectra steeper than around E~ 2 ^ above the 
vF v peak. Bursts such as the Superbowl one (GRB 930131) 
with high energy index (3 ~ —2 may be much more compatible 
with an SSC scenario. 

While this assessment may achieve closure in the 
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CGRO Composite Spectrum for GRB910503 
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CGRO Composite Spectrum for GRB910601 
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Fig. 4. — Photon spectra from the synchrotron self-Compton model for the bright bursts (a) GRB 910503 and (b) GRB 910601. The data compilations 
from the Compton Gamma-Ray Observatory (the operating energy bands of the three detecting instruments, BATSE, COMPTEL and EGRET, are again 
as indicated) are from Schaefer et al. (1998). The solid curves correspond to the resulting synchroton self-Compton spectrum that best approximates the 
observed continuum [labelled SSC 'fit']. Model parameters included 5 = 6.0, n = 2.0 and e ss 15 for GRB 910503, and 8 = 8.0, T) = 2.0 and e m 200 
for GRB 910601. In both panels the thermal (labelled T) and non-thermal (labelled NT) components are illustrated by the dotted and dashed curves, 
respectively. The small thermal contributions in both cases render the values of e quite uncertain, so that the depicted values should be regarded as 
approximate lower bounds to e. Adjustments along the energy axis yielded rf^,B±/B CT = 5.0 X 10 — 3 for GRB 910503 and F'y^B±/B CT = 3.3 X 10~ 3 
for GRB 910601, both implying sub-Gauss fields for typical values of r ~ 7 ~ 300 . 



Swift/GLAST era, a caveat is possibly provided if the syn- 
chrotron continuum is strongly self-absorbed. This is the sce- 
nario envisaged by Panaitescu & Meszaros (2000) and Liang, 
Boettcher & Kocevski (2003), with the latter group suggest- 
ing that SSC with moderate optical depths (r T ~ 0.1-0.4) 
could explain X-ray excesses seen in some burst spectra. With 
such strong absorption, the supply of soft photons for inverse 
Compton upscattering is more sharply peaked than even the 
synchrotron continuum from electrons possessing a truncated 
power-law distribution, thereby approximately mimicking the 
pure inverse Compton spectral shapes explored in Section 13.31 
The requirement that self-absorption be effective in narrowing 
the synchrotron self-Compton peak is that the turnover energy 
£t exceeds that of the synchrotron peak, which is below the 
X-ray band. Using Eq. J5J, one quickly arrives at the requisite 
condition: 

V7t & u{8) £ -?- , (26) 



Oi{ Speak 



for 



U (6) = 2-(^)/ 2 V3 (( 5-i)r(| + l)r(| + ^ 



(27) 



Here e peQ fe < 1 is the energy of the vF v peak in the prompt 
emission. For 2 < S <; 5 , one deduces that ui(5)/a{ ~ 110- 
160 , so that optical depths r T < 1 can be tolerated for typical 
r ~ 300 without pushing j T too low, i.e. into the mildly- 
relativistic regime. Detailed computations of self-absorbed SSC 
spectra are deferred to future work. 

3.6. Broadband Observational Constraints 

In addition to the soft gamma-ray spectral shape, a probing 
diagnostic on the emission mechanisms is whether their extrap- 
olation down to the optical band can be accommodated by op- 
tical monitoring data for bursts. Such broadband constraints 



turn out not to be profoundly limiting for synchrotron or pure 
inverse Compton emission. However, they do impact SSC mod- 
els, since these can potentially have a prominent synchrotron 
component in the optical that might violate upper bounds or, 
in the case of GRB 990123, prompt optical detections. It is 
straightforward to assess how optical band data can constrain 
SSC model parameters using the formalism derived here. 

First, the focus is on observational data. GRB 990123 (with 
redshift z — 1.6 ) is to date the only burst with a prompt optical 
detection (Akerlof et al. 1999), peaking at around magnitude 
8.8 some 47 seconds after BATSE trigger and contemporaneous 
with the BATSE signal. These prompt observations were made 
by the Robotic Optical Transient Search Experiment-I (ROTSE- 
I). The recent bright burst GRB 030329 (redshift z = 0.168) 
was only detected optically by ROTSE-III in the afterglow epoch 
at visual magnitude 12.5 over an hour after the burst trigger on 
HETE-2 (Smith et al. 2003). Kehoe et al. (2001) report up- 
per bounds of around magnitude 13-15 to half a dozen bursts, 
in two cases starting just 15-20 seconds after trigger. These 
low limiting magnitudes can potentially be severely constrain- 
ing on models with a high optical to gamma-ray flux expecta- 
tion, though contemporaneous multiwavelength detections are 
requisite for unambiguous diagnostics. 

Taking magnitude 9 as a benchmark optical signal for the 
purposes of exploring broadband constraints, a number that 
must be considered modulo the source distance and intrinsic 
luminosity, this corresponds to a differential photon flux of the 
order of 10 6 photons sec -1 cm -2 keV -1 at around 2eV. This 
can be compared with the BATSE spectrum of GRB 910503, 
a burst brighter than most of those in the sample of Kehoe et 
al. (2001), which has a time-integrated differential flux of 0.1 
photons sec -1 cm~ 2 keV -1 at around 100 keV. The interpola- 
tion between such an optical magnitude and gamma-ray flux is 
a roughly e -15 spectrum. Immediately, one infers that syn- 
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chrotron and inverse Compton models for the gamma-rays can 
accommodate such an optical magnitude; only limiting mag- 
nitudes of around 19 or higher can provide conflicts for these 
mechanisms. Note that a similar conclusion would apply for 
pitch angle synchrotron models, which possess very flat low en- 
ergy asymptotic spectra like inverse Compton scenarios. 

The situation for SSC broadband spectra is very different. 
The ratio of the energy of the SSC peak to that of the syn- 
chrotron peak is simply of the order 47^ in , which would set 
the synchrotron peak in the infra-red for Lorentz factors of 
7min ~ 10 3 , independent of the bulk flow T . Lower 7 m ; n 
can move the synchrotron peak into the optical and even the 
UV. The ratio of emissivities at the synchrotron and SSC 
peaks is quickly determined from Equations 1)12(1 and l(21|l to be 
N syn /N SS c ~ 87^ in /(r T /e) in the notations of this paper, i.e. 
with fg — ► e7 T 7/ 1 ~' s /(S — 1) . Hence, for r T fg ~ 1 , the peaks 
of the two components would lie roughly on an e _1 spectrum. 
Moreover, opting for 7 m ; n ~ 300 to position a synchrotron peak 
in the optical window, normalizing to the BATSE spectrum of 
GRB 910503 would yield a count rate an order of magnitude be- 
low the magnitude 9 benchmark posited above. At face value, 
this would render the SSC scenario compatible with limiting 
magnitudes of around 11 for this source. Yet Comptonization 
of the spectrum is not evident in the gamma-rays, so it is in- 
ferred that r T < 1, thereby increasing the ratio N syn /N ssc and 
lowering the permissible limiting magnitude for a viable SSC 
model. Accordingly, for the only source with a prompt optical 
detection, GRB 990123, for the SSC model not to overestimate 
the observed optical flux, the synchrotron peak has to move 
somewhat out of the optical, though not much. This amounts 
to a modest constraint in ( / y m i n , r T , 5) phase space that can 
be quickly deduced, though at this juncture is not particularly 
enlightening; such an exercise will be more timely once Swift 
is providing a number of simultaneous hard X-ray and optical 
detections of bursts. 



4. CONCLUSION 

The interpretative modeling presented here provides a strong 
indication that the synchrotron and inverse Compton processes 
may provide a more viable explanation of prompt burst spec- 
tra than the SSC mechanism. The spectroscopic bias against 



SSC is due to its inherently broad character around the vF v 
peak, and is only marginally constrained by broad-band consid- 
erations, i.e. limiting magnitudes in the optical band. However, 
strong self-absorption of the synchrotron component can provide 
an escape clause for SSC scenarios, though this imposes signifi- 
cant fine-tuning of parameter space. Inverse Compton provides 
an attractive possibility for the vast majority of bursts, as does 
pitch angle synchrotron emission, though quasi-isotropic syn- 
chrotron radiation may only be reconciled with around 2/3 of 
bursts at the lowest energies in the BATSE window. Discrim- 
ination between these processes may prove possible after Swift 
launches and in the GLAST era, though polarimetric measure- 
ments would enhance the diagnostic capability of observations. 
In analyzing SSC spectral models, attractively compact analytic 
forms for the synchrotron self-Compton cmissivity from thermal 
and truncated power-law electron distributions were developed; 
these will prove useful for applications of SSC spectral models 
to a variety of astrophysical environments. 

A distinctive conclusion of this work is that agreeable "fits" 
with either synchrotron or inverse Compton emission scenar- 
ios are only attainable if the underlying electron distribution is 
profoundly dominated by a non-thermal component. This result 
is in contradistinction with standard understanding of accelera- 
tion mechanisms such as the diffusive (Fermi) process at shocks, 
where the non-thermal particles are drawn probabilistically from 
the shocked thermal population, and therefore are subdominant 
in number. This constraint has profound implications for any 
proposed GRB model, providing a challenge for theorists in an 
era when computational resources are permitting new avenues 
for exploration using intensive numerical simulations of particle 
acceleration at shocks. 
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APPENDIX 



REDUCTION OF THE SYNCHROTRON SELF-COMPTON RATE FOR NON-THERMAL ELECTRONS 

In this Appendix, developments leading to the form in Eq. 11201 for the rate of synchrotron self-Compton emission from a truncated 
power-law distribution of electrons are expounded. The starting point is the insertion of i C sc times the synchrotron rate in Eq. (|11[> . 
serving as the soft photon density, directly into the expression in Eq. (|16(l for the inverse Compton emissivity. This presents a triple 
integral that can be reduced to a compact single integration by appropriate manipulations. The Lorentz factor 7 of the electron 
that upscatters the synchrotron photons proves to be a less convenient integration variable than z = e/(47 2 e s ) , which is defined in 
Eq. ijT7|l . Likewise, to proves expedient to use 

Q = *% B ". , , (Al) 

3-B y7min sin 9 

as an integration variable, rather than the synchrotron photon energy e s 
to an SSC emissivity of 



Performing these two changes of variables leads quickly 
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where N ssc is the SSC rate scale factor given in Eq. (JSIJ, and x is the dimensionless SSC energy variable defined in Eq. (J55J. 
Here, £ and the maximum value z max of z can be written 



£ = max{l, g} , z n 



min< 1, 






(A3) 



The z -integration over the inverse Compton scattering kernel leads to a definition of g$ : 

9t(z) = [ Z f(<l)q i5 - 1)/2 d q = zW*G s (z) , 
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The q integration can be reduced by noting that the definitions of £ and z max establish two critical values of q , namely q\ = x/lJ 
and q 2 — 1 . which delineate distinct dependences of the integrand on q . The relative ordering of q\ and q 2 depends on the value 
of y assumed: y > x =► 92 > Qi and y < x ^ 12 < qi ■ This naturally divides the y integration into two ranges, [0, x] > and 
[x, oo) . The developments are straightforward, and are facilitated using an hs function defined via 
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gsjq) 



« 
Az 
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The y > x ca se requires an integration by parts in one contributing q interval, and one quickly arrives at the compact form 
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i.e. Eq. i)20[l . using the following functional forms (with G5 and i/^ defined as just above): 
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for z->l/qi = y/x , and for z -> g x = x/y 
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Observe that as z — > 1 , P«5 — ► Q5 , as demanded by continuity of the spectrum. 

Asymptotic dependences of the SSC rate on the upscattered photon energy, encapsulated in the proportionalities of Eq. (|23|l . can 
be simply derived. In the limit x "C 1 , both the Pa and <5«5 terms contribute to Eq. I|A6J1 . with the integrals being dominated by 
the i/«l domain. The Bessel function can be replaced by its power-law asymptotic form K 5 / 3 (z) « 2 2 / 3 r(5/3) y~ 5 ^ 3 , resulting 
in elementary integrals. For 5 > 1/3 the integration converges to 
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The x S> 1 case is less simple, though the integration over Qs does not contribute significantly since it is exponentially suppressed 
by the dependence of the Bessel function. The P5 integration can be developed using identity 6.561.16 of Gradshteyn and Ryzhik 
(1980), which can be specialized to the case appropriate for this analysis: 
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The derivative of this with respect to [i can also be used to expedite the developments, and can be written in the form 

I x" log e x K 5/3 (x) dx = \ I(ji) U U + f\ + v/f - ±) + 2 log e 2 1 , (All) 
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where "4>(x) = d/dx[log e T(x)} is the Digamma function. Setting /i 
by infinity, one readily arrives at the result 



(S + l)/2 , since the upper limit to the P$ can be replaced 
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X > 1 



The logarithmic term will clearly dominate for extremely large values of x ■ 

THE SYNCHROTRON SELF-COMPTON RATE FOR THERMAL ELECTRONS 

The comparison of synchrotron self-Compton (SSC) models invoking shock acceleration with burst spectra requires also the 
SSC rate for relativistically hot thermal particles, which implies using a four-dimensional integral defined by inserting the thermal 
contribution to Eq. (JjjJ into Eq. J3J and then the result subsequently i nto E q. I|16|l . Rather than laboriously reducing the integral 
along the lines of the developments of Appendix A, the result in Eq. (|A6|) can be used to advantage. Taking the limit 5 — ► oo 
generates the SSC emissivity for monoenergetic electrons with Lorentz factor 7 m ; n . The result can then be convolved with the 
thermal electron distribution to yield the required emissivity. As S becomes very large, Pg(z) approaches zero at all positive z 
less than unity, so that the integration over Pg(z) in Eq. I|A6|I does not contribute. Using the result 



lim SQs(z) = j(z) = 2(z + 2)(z-l)-2(l + 2z) log e z 
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the integration over the thermal electron distribution produces a double integral that can be evaluated by reversing the order of 
integrations, yielding 
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The first integral defining J can be expressed as a sum of incomplete Gamma functions and one of their derivatives, using 
identities 3.381.3 and 4.358.1 of Gradshteyn and Ryzhik (1980), though this does not expedite the evaluation of J . However, the 
integral representation of the incomplete Gamma function in 8.353.3 of Gradshteyn and Ryzhik (1980) proves useful, leading to the 
alternative integral forms posited on the second line of Eq. I|B3|1 . These latter forms are readily amenable to numerical computation 
in a fashion similar to evaluation of Euler's integral form for the Gamma function. Note that the first form on the second line is 
expedient for k < 1 cases, while the second form on this line, which was obtained using the fact that the integrals of ve~ v s(v) 
and v 2 e~ v s(v) on [0,oo) are identically zero, is advantageous for evaluation when n ;> 1. Asymptotic forms for J are simply 

obtained: J{k) w 1184/ (225k) for k«1, and J{k) « 32e _K /K 4 for k>1. The appropriate asymptotic limits of Eq. 1B2J1 can 
then be quickly determined: 
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where the low argument power-law dependence for the Bessel function K 5 i 3 (y) is used when Xt *C 1 , and the method of steepest 
descents is employed to derive the Xt 3> 1 case, for which the exponential portion of K 5 / 3 (y) is sampled. 

Note that the SSC emissivity for a monoenergetic electron distribution can be determined by omitting the convolution with the 
thermal distribution, so that Eq. (|Bip can be used efficaciously to yield 
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As j(z) is a monotonically decreasing function of z on (0,1], for x ^ 1 it flattens the SSC spectrum relative to that for 
synchrotron emission, and so effects a broadening of the vF v peak. This monoenergetic formula can be convolved with arbitrary 
electron distributions to obtain resultant synchrotron emissivities. By the same token, Eq. l|B2|l can be used to facilitate a Laplace 
transform approach for special cases for the electron distribution. 



REFERENCES 



Akerlof, O, et al. 1999, Nature 398, 400. 

Band, D., et al. 1993, ApJ413, 281. 

Baring, M. G. 1988a, MNRAS235, 51. 

Baring, M. G. 1988b, MNRAS235, 79. 

Baring, M. G. 1990, MNRAS 244, 49. 

Baring, M. G. 1994, in Proc. of the Second Compton Symposium, eds. 

Gehrels, N. & Norris, J., (AIP Conf. Proc. 304, New York) p. 639. 
Baring, M. G. 1995, Proc. of the Erice NATO ASI Currents in High Energy 

Astrophysics, eds. Shapiro, M., Wefel, J. and Silberberg, R., (Kluwer 

Academic, Dordrecht) p. 21. 
Baring, M. G. 1999, in Proc. of the 26th International Cosmic Ray 

Conference, Vol. IV, p. 5, held in Salt Lake City, Utah. 
Baring, M. G. 2002, Pub. Astron. Soc. Aust. 19, 60. 
Baring, M. G., Ellison, D. C. & Jones, F. C. 1993, ApJ409, 327. 
Baring, M. G., Ogilvie, K. W., Ellison, D. C, & Forsyth, R. J. 1997, 

ApJ476, 889. 
Bckcfi, G. 1966, Radiation Processes in Plasmas, (Wiley and Sons, New 

York) 
Blandford, R. D. U Eichler, D. 1987, Phys. Rep. 154, 1. 
Brainerd, J. J. & Lamb, D. Q. 1987, ApJ313, 231. 
Brainerd, J. J. & Petrosian, V. 1987, ApJ 320, 703. 
Briggs, M. S., et al. 1999, ApJ 524, 82. 
Crider, A. & Liang, E. P. 1999, A&AS 138, 405. 
de Jager, O. C, et al. 1996, ApJ 457, 253. 

Dermer, C. D., Chiang, J. & Bottcher, M. 1999, ApJ 513, 656. 
Dermer, C. D. & Schlickeiser, R. 1993 ApJ 416, 209. 
Drury, L. O'C. 1983, Rep. Prog. Phys. 46, 973. 
Ellison, D. C, Baring, M. G. & Jones, F. C. 1996, ApJ 473, 1029. 
Ellison, D. C. k. Double, G. P. 2002, Astroparticle Phys. 18, 213. 
Ellison, D. C, Jones, F. C, & Reynolds, S. P. 1990, ApJ 360, 702. 
Ellison, D. C, Mobius, E., & Paschmann, G. 1990, ApJ 352, 376. 
Epstein, R. I. 1973, ApJ 183, 593. 
Epstein, R. I. & Petrosian, V. 1973, ApJ 183, 611. 
Forman, M. A., Jokipii, J. R. & Owens, A. J. 1974, ApJ 192, 535. 
Gallant, Y. A., Hoshino, M., Langdon, A. B., Arons, J. & Max, C. E. 1992, 

ApJ 391, 73. 
Ghisellini, G., Celotti, A. & Lazzati, D. 2000, MNRAS 313, LI. 
Giacalone, J., Burgess, D., Schwartz, S. J. & Ellison, D. C. 1992, 

Geophys. Res. Lett. 19, 433. 
Giacalone, J., Burgess, D., Schwartz, S. J. & Ellison, D. C. 1993, ApJ 402, 

550. 
Goodman, J. 1986, ApJ 308, L47. 
Gonzalez, M. M., et al. 2003, Nature 424, 749. 
Gosling, J. T., et al. 1981, J. Geophys. Res. 86, 547. 
Gradshteyn, I. S. and Ryzhik, I. M. 1980, Table of Integrals, Series and 

Products, (Academic Press, New York) 
Granot, J., Piran, T. & Sari, R. 2000, ApJ 534, L163. 



Hoshino, M., Arons, J., Gallant, Y. A. & Langdon, A. B. 1992, ApJ 390, 

454. 
Hoshino, M. & Shimada, N. 2002, ApJ 572, 880. 
Jackson, J. D. 1975, Classical Electrodynamics, (Wiley and Sons, New 

York) 
Jokipii, J. R., Kota, J., & Giacalone, J. 1993, Geophys. Res. Lett. 20, 1759. 
Jones, F. C. 1968, Phys. Rev. Lett. 167, 1159. 
Jones, F. C. & Ellison, D. C. 1991, Space Sci. Rev. 58, 259. 
Jones, F. C, Jokipii, J. R. & Baring, M. G. 1998, ApJ 509, 238. 
Kang, H. & Jones, T. W. 1997, ApJ 476, 875. 
Kehoe, R., et al. 2001, ApJ 554, L159. 

Kucharek, H., & Scholer, M. 1995, J. Geophys. Res. 100, 1745. 
Liang, E. P., Boettcher, M. & Kocevski, D. 2003, in Gamma-Ray Burst 

and Afterglow Astronomy 2001, Eds. G. R. Ricker & R. K. Vanderspek 

(AIP Conf. Proc. 662), p. 295. 
Liewer, P. C, Goldstein, B. E. & Omidi, N. 1993, J. Geophys. Res. 98, 

15,211. 
Lloyd, N. M. & Petrosian, V. 2000, ApJ 543, 722. 
Medvedev, M. V. 2000, ApJ 540, 704. 
Meszaros, P. 2001, Science 291, 79. 

Meszaros, P., Rees, M. J. & Papathanassiou, H. 1994, ApJ 432, 181. 
Pacynski, B. 1986, ApJ 308, L43. 
Panaitescu, A. & Meszaros, P. 2000, ApJ 544, L17. 
Pavlov, G. G. & Golenetskii, S. V. 1986, Ap&SS 128, 341. 
Piran, T. 1999, Phys. Rep. 314, 575. 
Preece, R. D., et al. 1998, ApJ 506, L23. 
Preece, R. D., et al. 2000, ApJS 126, 19. 
Preece, R. D., et al. 2002, ApJ 581, 1248. 
Rees, M. J., & Meszaros 1992, MNRAS 258, 41. 
Rybicki, G. B. & Lightman, A. P. 1979, Radiative Processes in 

Astrophysics, (Wiley and Sons, New York) 
Sari, R. & Esin, A. 2001, ApJ 548, 787. 
Schaefer, B. E., et al. 1998, ApJ 492, 696. 

Scholer, M., Trattner, K. J. & Kucharek, H. 1992, ApJ 395, 675. 
Shimada, N. & Hoshino, M. 2000, ApJ 543, L67. 
Shimada, N., et al. 1999, Ap&SS264, 481. 
Smith, D. A. et al. 2003, ApJ 596, L151. 
Tavani, M. 1996a, Phys. Rev. Lett. 76, 3478. 
Tavani, M. 1996b, ApJ 466, 768. 
Terasawa, T., et al. 2001, in Proc. 27th International Cosmic Ray 

Conference. (Copernicus Gesellschaft. Germany), SH2.3/SH119 [see 

http://www. Copernicus, org /icrc/ program, htm 
Trattner, rv. J. & Scnoler, M. lyyi, ueopnys. ties. Lett. 18, 1817. 
Trattner, K. J. & Scholer, M. 1993, Ann. Geophysicae 11, 774. 
Westfold, K. C. 1959, ApJ 130, 241. 
Zhang, B. & Meszaros, P. 2001, ApJ 559, 110. 



